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Abstract — In this paper, a joint power allocation algorithm with 
minimum mean-squared error (MMSE) receiver for a cooperative 
Multiple-Input and Multiple-Output (MIMO) network which 
employs multiple relays and a Decode-and-Forward (DF) strategy 
is proposed. A Distributed Space-Time Coding (DSTC) scheme 
is applied in each relay node. We present a joint constrained 
optimization algorithm to determine the power allocation pa- 
rameters and the MMSE receive filter parameter vectors for 
each transmitted symbol in each link, as well as the channel 
coefficients matrix. A Stochastic Gradient (SG) algorithm is 
derived for the calculation of the joint optimization in order to 
release the receiver from the massive calculation complexity for 
the MMSE receive filter and power allocation parameters. The 
simulation results indicate that the proposed algorithm obtains 
gains compared to the equal power allocation system. 
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I. Introduction 

MIMO wireless communication systems employ multiple 
collocated antennas in both source and destination node in 
order to obtain the diversity gain and combat multi-path 
fading. The different methods of STC schemes, which can 
provide a higher diversity gain and coding gain compared 
to un-coded systems, are also utilized in MIMO wireless 
systems. However, it is often impractical to apply MIMO 
in mobile communication systems due to the high cost and 
the size of mobile terminals. Cooperative MIMO systems, 
which employ multiple relay nodes between the source and 
destination node as the antenna array, apply distributed di- 
versity gain in wireless communication systems Q. Among 
the links between the relay nodes and destination nodes, 
cooperation strategies, such as Amplify-and-Forward (AF), 
Decode-and-Forward (DF), and Compress-and-Forward (CF) 
E) and various DSTC schemes in 0, O and can be 
employed. 

Recent contributions in the cooperative communications 
area lie in the power control problem in transmitters using the 
AF strategy lf7l- lfTTI . A central node which controls the trans- 
mission power for each link is employed in [7 |. Although the 
central control power allocation can improve the performance 
significantly, the complexity of the calculation increases with 
the size of the system. The works on the power allocation 
problem for the DF strategy measuring the outage probability 
in each relay node with single antenna and determining the 
power for each link between the relay nodes and destination 
node, have been reported in [4|-[6|. The diversity gain is 
sacrificed by the utilization of a single antenna in relay nodes. 



In this paper, we propose an adaptive power allocation 
algorithm with linear MMSE receiver for cooperative MIMO 
systems employing multiple relay nodes with multiple anten- 
nas to achieve a DF cooperating strategy. The power allocation 
parameters and the linear MMSE receive filter parameter vec- 
tors can be determined and fad back to each transmission node 
through a feedback channel that is error free and delay free. 
The joint estimation algorithm for power allocation parameters 
in each link and the MMSE receive filter for each symbol 
is derived. By utilization of an SG algorithm from [1|, the 
complexity of the calculation will be decreased compared with 
the MMSE-based expressions that require matrix inversions. 
Also the channel estimation is done by using the SG algorithm 
before the computation of the power allocation parameters. 

The paper is organized as follows. Section II provides the 
multi-hop cooperative MIMO system with multiple relays 
applying the DF strategy and DSTC scheme. Section III 
describes the constrained power allocation problem and linear 
MMSE detection method, and in Section IV, the proposed it- 
erative SG algorithm is derived. Section V focus on the results 
of the simulations and Section VI leads to the conclusion. 

II. Cooperative System Model 

The communication system under consideration, 5 shown 
in Fig. 1, is a MIMO communication system transmitting 
through a multipath channel from the source node to the 
destination node. There are n r relay nodes, applying a Decode- 
and-Forward (DF) scheme as well as space-time coding (STC), 
between the source and the destination node, and N antennas 
at each node for transmitting and receiving. A multi-hop com- 
munication system can be achieved by broadcasting symbols 



from the source to n r — th relay nodes as well as to the 
destination node in the first phase, followed by transmitting 
the detected and re-encoded symbols from each relay node 
to the destination node in the other phases. We consider only 
one user at the source node in our system that has N Spatial 
Multiplexing (SM)-organized data symbols packed in each 
packet. The received symbols at each relay and the destination 
node, denoted as rsR„ T and r$D respectively, are detected 
by a linear MMSE receive filter at each receiving node. We 
assume that the synchronization in each node is perfect. The 
transmission between the source node and the k — th relay 
node, and the destination node can be described as follows 

tsrM = F k [i\AsR k [i\s[i\ + n S R„[i], (1) 
r SD [i] = H[i]A SD [i\s[i] + n SD [i], (2) 
i = l,2, ... ,N, k = 1, 2, ... n r 

where the N x 1 vector nsR k [i] and ns£>[i] denote the zero 
mean complex circular symmetric Additive White Gaussian 
Noise (AWGN) vector generated in each relay and the desti- 
nation node with the variance of a 2 . The transmitted symbol 
vector s[i] = [si[i], s 2 [i], sn^]] contains N parameters, 
and has a co variance matrix = a 2 1, where E[-] 

stands for expected value, (-) H denotes the Hermitian operator, 
a 2 is the signal power which we assume to be equal to 1 and J 
is the identity matrix. Fk[i] and H[i] are the N x N channel 
gain matrix between the source node and the k — th relay 
node, and between the source node and the destination node, 
respectively. .As£>[z] and Asij fc [i] are the diagonal N x N 
power allocation matrices with complex parameters asc[j] 
and asR k [i] assigned to each symbol vector s[i] transmitted 
to the k — th relay node and the destination node. 

After filtered in each relay node, the detected symbols will 
be re-encoded by a N x T distributed space-time coding 
(DSTC) matrix and assigned a power allocation parameter 
matrix and then, forwarded to the destination node. Notice that 
only the relays which can detect the received symbols correctly 
will forward the encoded symbols to the destination node 
because the interference between the received symbols after 
space-time coding will be increased if the encoded symbol 
vectors are different. Define n rc ii a bi with length L to be the 
relay set which can implement the correct detection. Then the 
relation between the I — th relay and the destination node can 
be described as 

T 

R Ri D[i\ = ^2G RlD [i]A t [i]m Rl Dt [i] + N Ri d W j 0) 
t=i 

/ £ ^reliable 

where the N x 1 matrix mjj,D t [*] is the t — th column of 
the DSTC matrix, I is the number of reliable relay nodes 
which can detect the received symbols correctly, and A t [i] is 
the diagonal matrix contains the power allocation parameter 
assigned to the t — th column of the re-encoded matrix. The 
N x N channel gain matrix is denoted by Gr^Ii], and the 
N xT AWGN matrix Nr^Ii] is generated in the destination 
node with variance a 2 . The N x T received symbol matrix 



Rrid[i\ in (3) can be transformed and expressed as a NT x 1 
vector rR 4 D[i] given by 

N 

[i] +n RlD [i], (4) 

where the T x 1 vector a jRiD [i] = [ai RjD [i], a TRlD [i]} is 
the power allocation parameter vector assigned for the j — th 
symbol Sj [i] . The diagonal NT x N matrix D j R D [i] = 
dia,g[dj 1 [i],dj 2 [i], ...,dj T [i]] stands for the effective channel 
coefficient matrix combined with the DSTC scheme and the 
channel matrix Gr^Ii]- After rewriting r^ ; £>[z] we can 
consider the received symbol vector in the destination node 
as L + 1 parts, one is from the source node and the remaining 
L are from the reliable relays, and write the received symbol 
vector for cooperative detection as 

J2f=i h 3 [i} a 3SD [i] + n SD [i] 
E^i D j RlD [*K' RlD [i] + n Rl D[i] 

. Ej=i d Jr 2 d \i]aj R2D [i\s jR2D [t] + n R2D \i] 



_ Y,j=i D 3R LD [i]aj RLD \i\sj RLD [i}+nR LD [i] _ 

N 

= [*]<**•[*] «j [*]+»»[*] 

3 = 1 

(5) 

where h 3 \i\ is the j — th column of the N x 
N channel coefficients matrix H[i] and Sj R[D [i] = 

™f SRl (Ef=i fj SRl \ihjsR, + n S R,[t\) is the detected 

symbol in the I — th relay node, where /,•__ [i] is the 
j — th column of the N x N channel coefficients matrix 

FsrSY The ( LT + l ) N x ( LN + 1) diagonal matrix 
Bj[i] = diag[h j,D jRiD \i], ... , D jRhD [i]] contains the 
channel gain elements of all the links between the reliable 
relays and the destination node. The (LN + 1) x 1 power 
allocation parameter vector a 3 \i\ = [aj sn [i], a,j R D [i]] T = 
[ a js[i],a jRlD \i], ... , a JRLD [i]] T , where a 0RD [i] denotes 
the power assigned for the DSTC matrix for I G n r ciiabic, 
containing all the power elements in each link. 

III. Joint linear MMSE receiver design with 
Power Allocation 

The MMSE receiver design with power allocation for 
every link between the source and the destination node 
as well as the reliable relays is derived as follows. If 
we define a (LT + 1)N x 1 parameter matrix Wj[i] = 
[wiji^tuiji], ...,w 1l+1 [i\] for i = 1,2, ...,N K /N to deter- 
mine the j - th symbol [i] and a«j R|D [i] for j = 1,2, ...,N, 
by using (5), the MMSE problem with power allocation can 
be described as 

[wj t0pt [i\,a jt0pt [i\} = arg min E[\\ Sj [i] - wf \i}r\i}\\ 2 } 

■w 3 [t], aj [i} 

(6) 



subject to 

JV n r N 

J2 \i]°<*j SD + a lsn k [i]<XjsR„ W = Pt ' 

j=l fe=lj=l 
L N 

i=i j=i 

where Pt is the total transmit power for the source node and 
all the reliable relay nodes. The joint linear MMSE receiver 
design and power allocation problem adjusts the receive filter 
w j,opt[i] an d the power allocation parameters aj, opt [i] with a 
power allocation constraint. It can be transformed into an un- 
constrained optimization problem using Lagrange multipliers 
[1] which has the following lagrangian 

L=E[\\ Sj [q-wfmm 2 }+xi(f2^sA^* jSD [i\ 

n r N 

+EE a W*kW*]- PT )+ (7) 

fc=l j=l 
L N 

J2 W a j Rl o [*] - p t), 
i=i j=i 

From (7), we can see there are three power allocation param- 
eters to be determined, which are power assigned to the link 
between the source and the destination node Q.j SD [i], between 
the source and each relay node ctj SR [*] an d between the l—th 
reliable relays and the destination node ctj R D [i]. By fixing 
these Qj [i] and taking gradient terms of (7) and equating them 
to zero we obtain the expression of w^ opt [i\ which is 

■Wj,o P t[i\ = R^PwAi]' W 

where R Wj ^ is the (LT + 1)N x (LT + 1)N covariance 
matrix equals to E[r 3 ;[i]rf[i\\ and the (LT + 1)N x 1 cross- 
correlation vector p w m = E \r 3 [i\sj . In the expression of 
R Wj [i] and p w ^ there exists the power allocation parameters. 
The expression of power allocation parameters can be obtained 
by fixing Wj. opt [i] in (7) and taking gradient terms with respect 
t0 a isD,opt[i], a jRiDt0pt [i\ and a 3SR ^ opt [i\ and equating them 
to zero, 



a 3SD,opt [i] 






(9) 


a 3R l D,Opt [i] 


= RJ 




(10) 


a jsR k ,opt [i] 


= R 1 

a iSR k 


[»] Pa *SR k I*]' 


(11) 



where the covariance matrices are 

R ajSD[i] = E [w? [t] hj [i] s 3 - [i] s* [i] hj [i] H w h [t]] + Ai 

R oiR lD W = E i w fi i { \ d 3R l o\i}sj[i\s*[i\df R[D \i}w jl [»]] + A 2 

R "i S R k H = E [ w fi ['K^is^WS^K^ WisR k ) s M 
(D jRlD [^ jRiD (wf SRk \i}f jsRk )s J [i]) H w Jl [i]] 
+ Ai 



and the cross-correlation vectors can be calculated as 

p ^ ;D w =^[< D M^H^W^W] 

^ SRfc H =^[« [i]£>w lD [i]aM I D 

(«i« fc W/j« t ))^[»>i[»1] 

where lOj, [i] is the MMSE receive filter parameter vector for 
the I — th receiver symbol vector n [i] . 

The expression of the MMSE detection vector and the 
power allocation elements depend on each other as well as 
the effective channel matrix Dj[i], and should be determined 
by iterating with initial values and the Lagrange multiplier Ai 
and A2 to obtain the result. 

IV. Adaptive Estimation Algorithm for MMSE 
Design with Power Allocation 

The formulas in the previous section describe the method 
to calculate the MMSE detection vector Wj [i] and the power 
allocation parameter otj[i] for each transmit symbol Sj[i], 
which require matrix inversions with high complexity as well 
as channel estimation in the iteration calculation. In this 
section, an adaptive estimation algorithm based on an SG 
algorithm will be presented to determine the MMSE receive 
filter, the power allocation parameters and the effective channel 
matrix without the inversion calculation. 

A. Adaptive SG Estimation for MMSE Receive Filter and 
Power Allocation 

In this subsection, we will present the adaptive SG esti- 
mation algorithm for the MMSE receive filter w 3 \i\ and the 
power allocation parameter a. 3 \i\. The problem is described in 
formula (6) and by using the Lagrange multiplier method [1] 
we can obtain the expression (7) which indicates the parame- 
ters depend on each other. As a result, we will develop a SG 
joint estimation algorithm with low complexity calculation to 
solve the problem in (7). 

Considering the Lagrangian function in (7) and computing 
the instantaneous gradient terms of it with respect to Wj [i] and 
a.j\i\, respectively, lead us to the equations (12)-(15), where 
dj[i] is the j — th column of the effective channel matrix 
Dj Ri D [i] with dimension TN x 1, and ej [i] is the error symbol, 
which indicates the distance between the transmitted symbol 
and the detected symbol, calculated by Sj[i] — to^[i]r[z] = 

Sj[i] - wf[i](^2j =1 Bj[i]a,j[i]sj[i] + n[i]), and (•)* denotes 
the conjugation. Notice that we need to determine the MMSE 
receive filter w j h sr [i] at each relay node in order to determine 
if the k — th relay can detect the received symbols correctly 
and then forward it to the destination node. We can calculate 
the instantaneous gradient terms of L with respect to Wj SRk [i] 
to obtain 

Vi TO » w - - w% [i\Dj [i)a 3R [i](F k [i]A SRk [i\s[i] 

JbH k 

N 

+ n S R k \i}) (sj [i] - wf [i] B 3 [i] a 3 - [i] s 3 [i] 
+ n[i]))* 

= - wf [i\Dj [i]a 3R [i]r SRk [i]e* [i], 

(16) 



3sn 1 ' 



Vi„,;[i] = -(Y. B M a M s iV\ +n[i\)(sj[i\ -wf [i](^2B 3 [i]a 3 [{\s 3 [{\ +n[i]))* = -r\i]e*\i], 



3=1 



3=1 



AT 



2 = 1 
JV 



2 = 1 



(12) 



-s*[«]hf [i](sj[i] - wf [i](^ ^ W^W^ W + «[*])) + a i«kd [»] = -^W^f W^ji WejW + A i a js D W ( 13 ) 



-s*[i]df [t\( Sj [t\ - wf B,-[i] 0j + n[i\)) + \2<* jRiD \i] = - a J[»]df [i]e,[i] + A 2 cy SlD [i] (14) 



R k (15) 

= Vt D jR,D WojujD W/i ss . + A2a iS|D M 



We can devise an adaptive SG estimation algorithm by using 
the instantaneous gradient terms of the Lagrangian which were 
previously derived with the SG descent rules [1] as 

(17) 

[il (18) 



N 



Wj[i + 1] 



Wi i 



w 



jsR k [* + !]= w j S R k [i] - A*V£„ 
aj SD [* + !] = a JSD [«] - 7V£ a ; 

<W [» + 1] = «J fl , D H - 7Vi Q ; 



,[< + !] 



[*] - 7Vi Q - 



(19) 
(20) 
(21) 



, where fi and 7 are the step sizes of the recursions for the 
estimation of MMSE parameter vectors which have to be de- 
termined before the estimation. The complexity of calculating 
Wj[i] is (O(N)) and (O(N)) for calculating aj SD [i} and 
a jR,n[i]' an( l (0(TN)) for aj SR [i], which are much less 
than that of the algorithm we described in Section III. As 
mentioned in Section II, all the MMSE receiver filters and 
power allocation matrices will be transmitted back to the relay 
nodes via a feedback channel which is assumed to be error- 
free in the simulation; however, in practical situation the errors 
at each relay node should be considered due to the property of 
broadcasting and the diversification of the feedback channels 
with time changes. 

B. Adaptive SG Channel Estimation 

In this subsection we will derive an adaptive SG algorithm 
for estimating the effective channel matrix B j [i] . The channel 
estimation can be described as an optimization problem 



are min E\\\ 
BAA L 



Bi 



i d, J Si 1 



(22) 



Define the received symbol vector r[i] = [ri [i], r^+i [i]]. 
The optimization problem described in (22) can be divided 
into three parts, which correspond to individually computing 
H[i], FsR k [i] and Dj R D [i]. We then derive an SG algorithm 
by calculating the cost function Ch, Cf and Co 

|21 



C H = E[\\r 1 \i}-H\i}As D \i]s[i\ 



(23) 



N 



C F =E[\\n\i] -J2D jRiD \i]a jRiD \i]wf SRi MFsrA 
A SRl [^+n SRl [i})\\ 2 ] 



(24) 



C D = E[\\n\i] - J2 D m,D [i\a, jRlD \i]s jsRi [i] 



(25) 



and then by taking instantaneous gradient terms of Ch, Cf 
and Cd with respect to H[i], F[i] and the jth column of 
effective channel vector Dj R D [i] 

VC H . {i] = -AsD[i\s\i}( ri \i] - H\i}AsD\i}s\i}), (26) 

JV 

V °f* M = - E D m lD [^a jRiD [i\wf SRi [i\(A SRl [»]«[»]) 



N 



(ri [*] ~ D iR l o\i]a jRlD \i]wf SRi \i]r SRl ), 



3=1 



N 



(27) 



VC D .[i] = -a jRiD [i]sj[i](ri[i] -^D jRiD [i]a jRiD [i] 

2 = 1 



°JSR 



,M)> 



(28) 

Considering the SG descent rules [1] and the result of gradient 
terms we can obtain the adaptive SG channel estimation 
expression which is 



H[i + 1] = H[{\ - 0VCH.fr 
F[i + l] = F[t\-F7C F .u i 



(29) 
(30) 



D jRiD [i + 1] = D jRiD [i] - WC D , RiD[i] (31) 

where (3 is the step size of the recursion. The adaptive SG 
algorithm for effective channel H[i], F[i] and Dj RD [i] re- 
quires the calculation complexity of (0(N 2 )), (0(TN)) and 
(0(TN)) and can determine the channel matrix accurately. 

V. Simulations 

The simulation results are provided in this section to assess 
the proposed algorithm. The system we considered is a DF 
cooperative MIMO system with different distributed STC 
schemes using QPSK modulation in quasi-static block fading 
channel with Additive White Gaussian Noise (AWGN), as 
derived in Section II. The bit error ratio (BER) performance of 
the joint power allocation using a linear MMSE receive filter 
(JPA-LMF) algorithm and the equal power allocation using 
a linear MMSE receive filter (EPA-LMF) algorithm with the 
power constraint employs different number of relay nodes and 



different STC schemes are compared. In the simulation we 
define the power constraint Pt as equal to 1, and the noise 
variance a 2 for each link is equal to 1 as well. 

The proposed JPA-LMF algorithm is compared with 
the EPA-LMF algorithm using the distributed-Alamouti (D- 
Alamouti) STBC scheme in [14] with n r = 1,2 relay nodes 
in Fig. 2. The number of antennas N = 2 at each node and the 
cooperative decoding delay at the destination node is T = 2 
time slots. The results illustrate that the performance of EPA- 
LMF algorithm is close to that of JPA-LMF algorithm when 
using the same DSTC scheme in lower E^/Nq circumstance; 
however, with the E^/Nq increase, the JPA-LMF algorithm 
obtains about 5dB of gains compared to the EPA-LMF algo- 
rithm using the same DSTC scheme to achieve the identical 
BER. The performance improvement of the proposed JPA- 
LMF algorithm is achieved with more relays employed in the 
system as an increased spatial diversity is provided by the 
relays. Fig. 3 illustrates the performance of the proposed JPA- 
LMF algorithm using different STC schemes under the same 
Eh /No condition and which employs n r = 1 relay node. The 
cooperative detection is achieved using a linear MMSE receive 
filter at the destination node. The randomised-Alamouti (R- 
Alamouti) STBC |fl5l is employed in the simulation and the 
number of antennas is set to N — 2 in each node. The 
simulation results shown in Fig. 3 indicate that the proposed 
JPA-LMF algorithm using R-Alamouti STBC obtains about 
5dB of gain compared to the EPA-LMF algorithm using 
the same DSTC scheme. There exists 3.5dB of gain of the 
simulation result for the R-Alamouti STBC utilizing the JPA- 
LMF algorithm compared to that for the D-Alamouti STBC 
utilizing the same algorithm to achieve an identical BER as 
illustrated in Fig. 3. 

VI. Conclusion 

We have proposed a joint power allocation and receiver 
design algorithm using a linear MMSE receive filter with 
the power constraint between the source node and the relay 
nodes, and relay nodes and the destination node. A joint iter- 
ative estimation algorithm for computing the power allocation 
parameter vector and linear MMSE receive filter has been 
derived. The simulation results illustrate the advantage of the 
proposed power allocation algorithm by comparing it with the 
equal power algorithm. The proposed algorithm can be utilized 
with different distributed STC schemes using the DF strategy, 
can be extended to the AF cooperation protocols and non- 
linear receivers lfj"8l . 
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